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A previously introduced real space renormalization-group treatment of the random transverse-field 
Ising spin chain is extended to provide detailed information on the distribution of the energy gap 
and the end-to-end correlation function for long chains with free boundary conditions. Numerical 
data, using the mapping of the problem to free fermions, are found to be in good agreement with 
the analytic finite size scaling predictions. 
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I. INTRODUCTION 

It has been known for some time that the properties 
of many random systems are dominated by rare regions. 
In random quantum systems at low or zero temperature, 
this effect is particularly pronounced. For example, in 
random quantum magnets, typical spin correlations at 
zero temperature can decay rapidly with distance while 
average correlations - which are measured by scattering 
experiments - decay much more slowly. In extreme cases, 
the average correlations and related quantities such as 
susceptibilities are dominated by rare pairs of spins whose 
correlations can be large even when they are very far 
apart. 

In the last few years, it has become clear that the sim- 
plest of all random quantum systems that undergoes a 
phase transition, the random transverse field Ising spin 
chain, exhibits these phenomena in a particularly dra- 
matic manner. Surprisingly, a tremendous amount of 
analytic information can be obtained about this system 
in the low energy, long length scale regime near the crit- 
ical point from a real space renormalization gppup (RG) 
analysis that yields exact asymptotic resultsEJ. So far, 
most of the results have concerned the behavior of av- 
erage quantities, such as the mean spin-spin correlation 
function and the susceptibility. But the general struc- 
ture of the RG implies that typical correlations are very 
broadly distributed and take an entirely differeni, form 
than the average correlations. Numerical studieso have 
yielded the distributions of spin correlations and energy 
gaps and found scaling behavior consistent with the gen- 
eral expected structure but there have, so far, not been 
analytic predictions for the distribution functions. 

In this paper, we focus on the behavior of random 
transverse field Ising chains with free boundary condi- 
tions and study the scaling behavior of the distributions 
of the energy gap and of the corcplations between the 
spins at opposite ends of the chainEl. We obtain analytic 
results for these distributions in the long-chain scaling 
limit and compare them with numerical results on finite 
length chains. The connections between the couplings in 



each particular sample and its properties are also studied. 

In addition to its own intrinsic interest, it is hoped that 
better understanding of correlations in this seemingly 
simple system may lead to progress on other more compli- 
cated systems - particularly higher dimensional random 
quantum systems with phase transitions which are not 
governed by RG fixed points accessible by conventional 
field theoretic methods. 



II. THE MODEL 

The model that we study is the one-dimensional ran- 
dom transverse-field Ising chain with Hamiltonian 



L-l 



(1) 



Here the {erf} are Pauli spin matrices, and the interac- 
tions Ji and transverse fields hi are both independent 
random variables, with distributions 7r(J) and p{h) re- 
spectively. Since the model is in one dimension, we can 
perform a gauge transformation to make all the Ji and 
hi positive. The numerical work used the following rect- 
angular distributions: 



p{h) 



1 for < J < 1, 
otherwise. 



h^^ 



'.Q for < /i < /lo, 
otherwise. 



(2) 



which are characterized by a single control parameter, /iq. 
The lattice size is L, and, in this paper, we will consider 
free boundary conditions. 
Defining 



A,i = In ft, 
A,7 ==1^ 



(3) 



where the overbar denotes-ian average over disorder, the 
critical point occurs whenQEl 



1 



A, 



(4) 



Qj = max{/ii, Ji}. 



(10) 



Clearly this is satisfied if the distributions of bonds and 
fields are equal, and the criticality of the model then 
follows from dualityt. 

We define a quantity, 6, which measures the deviation 
from criticality by 



6 = 



var h + var J 
where ly is given by 

Id 



var h + var J ' 



(5) 



(6) 



with "var" the variance over the distribution of the cou- 
plings in the Hamiltonian. 

For the distribution of Eq. (||), we see that from 
Eq. (^, the critical point corresponds to /iq = 1, and, 
from Eq. (|5|), (5 is given by 



while Eq. (g) gives 



(5 = - Inho, 



1. 



(7) 



(8) 



We will focus on the energy gap, AE, between the 
ground state and first excited state, and the end-to-end 
correlation function 



(9) 



This latter quantity, for reasons that are somewhat dif- 
ferent than those in the pure IsingiJiiodel, is rather easier 
to analyze than bulk correlationsB. We first summarize 
the RG method and the basic picture of the system that 
emerges, then present the analytic results, and finally 
present the numerical results and compare them with the 
analytic predictions. In the Appendix, the distribution of 
the energy gap with periodic (rather than free) boundary 
conditions is computed. 



III. ANALYTICAL RESULTS 

In this section we use the methods of Ref. |l| to obtain 
the scaling forms of the distribution of In AE and In Ci,l 
which will be valid in the limit L, — InAE and — InCi.L 
large and the distance to the critical point small. We 
will follow the notation of Ref. |^ and refer to equations 
in that paper by, e.g. F(5.28). 



A. Renormalization Group Analysis 

The basic renormalization group (RG) procedure from 
which the low energy, long distance behavior can be ob- 
tained is as follows: First choose the largest coupling in 
the system 



If the coupling is a field, say ft,2, then decimate the asso- 
ciated spin, leaving, by second order perturbation theory, 
an effective bond 



J 



J12J23 



13 



< /l2 



(11) 



which couples together the spins on the two sides of the 
decimated one. For our purposes here, it is important to 
keep track of the correlations, so we note that first order 
perturbation theory yields a correlation of the decimated 
spin with each of its neighbors, e.g.: 



J 



12 



H2' 



(12) 



If instead the strongest coupling is a bond, say J12, 
then join together the spins on both sides of this bond to 
make a spin cluster, and assign it an effective spin cluster 
variable ai2. The couplings between this cluster and its 
neighbors are the same as before; e.g. — Joi(Tf20'oi ^'^d 
the new effective transverse field acting to flip the spin 
cluster is 



112 



hih2 
1^ 



< J 



12 



(13) 



by second order perturbation theory. Since the two spins 
in the cluster will tend to flip together, their correlations 
will be strong 



J12 



J12 



(14) 



By decimating either a spin or a bond we have de- 
creased the number of degrees of freedom in the system. 
The resulting effective Hamiltonian has exactly the same 
form as before but we now have spin clusters and effective 
bonds with lengths {Ici} and {Isi}, respectively, which 
now vary. This decimation procedure is iterated with the 
largest of the couplings in the new effective Hamiltonians 
il = max{ft,i, Ji} gradually being reduced. 

This renormalization group transformation has two es- 
sential features. First, the effective bonds and fields re- 
main independent (although, e.g. the length of a bond 
and its strength are correlated). Second, if the system 
is near critical, the distributions of the basic logarithmic 
couplings: 



and 



C, = HQ/ J,) 



13, = 



(15) 



become broader and broader — even on this logarithmic 
scale — due to the additive nature of the recursion rela- 
tions for C and /3; e.g. from Eq ([III), 



2 



Cl3 — Cl2 + C2 



(16) 



This implies that at low energy scales, 17 ^ fi/, corre- 
sponding to 



Inn 



(17) 



large, the perturbation expansion used in deriving 
Eqs. (|ll-|l4) becomes very good as the neighboring cou- 
plings of a decimated coupling (which is by prescription 
equal to fi) will typically have magnitude <C ri. The oc- 
casional bad case in which two neighboring couplings are 
comparable can be shown (see F Section VI) to lead only 
to corrections to scaling and to non-universal amplitudes. 

The natural quantities to work with are the joint distri- 
butions of the logarithmic couplings and lengths P{C, Ic) 
and R{(3,Ib) of the bonds and clusters, respectively. 
The non-linear integro-differential recursion relations for 
these (F2.13) and F(2.14) are simply obtained. In terms 
of the Laplace transform variables. 



PiCy) 



die-y^p{C,i) 



(18) 



the recursion relations for each y decouple due to the 
simple additive properties of lengths. As shown in F Sec 
II, there is a special four parameter family (for each y) 
of exact solutions of the RG recursion relations. General 
well behaved initial distributions converge very rapidly to 
this special family and then at low energy scales converge 
to a subfamily of scaling solutions parametrized simply 
by 5, defined in Eq. (|^), and and a basic length scale, 
It], defined in Eq. (H). Note that the particular initial 
distributions we have used for our numerical studies are 
in the special family at the critical point but not away 
from it; this will result only in corrections to scaling. 



B. Obtaining End— to— End Correlations and Gaps. 

In order to obtain the properties of a finite chain of 
length L with free boundary conditions, we must follow 
the correlations of the first and last spins erf and (t£ with 
the active — i.e. undecimated — spins that remain at log- 
energy scale r. Since we will always be interested in the 
tj^ correlations, we will drop the z superscript for conve- 
nience. As discussed above and in detail in F Sec VB, 
the correlations of an end spin, say the first, cti, can be 
obtained by first observing that just after it is decimated, 
say at a scale we denote Fi, its correlations with the spins 
in the leftmost remaining cluster, which has cluster spin 
(T2(ri), are given by Eq ( p^ ) with the appropriate effec- 
tive bond at that scale J12 = e"^"^"^"^^. and effective 
field hi — e^^i on the just-decimated cluster which the 
(Ti was in, i.e. 



(19) 



^23 ^3 



FIG. 1. Schematic of renormalization group decimation 
and its use for obtaining end-to-end correlations. Clusters 
of spins (vertical lines) are demarcated by closed curves, deci- 
mated spins are indicated by small x's and decimated clusters 
by large X's. The dotted lines denote the effective bonds cou- 
pling the spin clusters. At a scale Fi, the cluster containing 
the first spin o"i was decimated with an effective bond strength 
e~'"i~''i^ coupling it to the left-most remaining cluster with 
spin (T2- Later, at scale r2 > Fi, the ct2 cluster was decimated, 
and the left-most remaining cluster is 0-3 as shown. The last 
spin, ol, is in a cluster that has not yet been decimated. 



Af(Fi) « C12. 



(20) 



Note that we have ignored here the corrections at earlier 
stages of the renormalization from Eq (|l^ as these will 
only give rise to multiplicative factors of order unity, cor- 
responding to additive factors of order unity in A, that 
arise from high energy scales Vl Vlj. These are only 
corrections to scaling. At a later stage of the RG, just 
after the leftmost cluster is itself decimated at scale F2 , 
its correlation with the new leftmost remaining cluster 
o'3(F2) is determined by the effective bond connecting 
them which has magnitude given by C23 ■ Therefore from 
Eq (0) 



Af(F2) = - ln(aia3(F2)) « Kf{Ti) + C: 



23 



(21) 



This process is illustrated schematically in Fig. ^ 

We can continue the process from both ends until there 
is only one remaining spin cluster with effective field 



hE -= e-^^ 



(22) 



The log-correlations of the first and last spins with this 
last cluster are given by the A's at this scale, Api^E) and 
Al{^e) (for "first" and "last") respectively. We thus 
have, for the end-to-end correlation function of interest. 



A = - In(Ci^L) = AFiTE) + AdTE) 



(23) 



with additive errors typically of order unity. The energy 
gap is simply given by the scale F^; at which the last 
degree of freedom — this last cluster — is decimated: 



G = - ln(A£;) « Ti 



(24) 



with 



In order to obtain the distributions of these quantities 
of interest, we must keep track of the joint distribution at 
scale F of the number of remaining clusters, N, their N 
log-effective fields {A}, the A^ — 1 log-bond strengths {d} 
and the contributions to the log-correlations of this stage 
of the RG, Af(F) and Al(F) from Eq (^), etc. A natural 
first guess is that, as in the infinite system, each of the 



3 



couplings will be independent. But, as shown explicitly in 
F Section IID, this is not correct: correlations arise from 
the constraint that the total length is L. For example 
an anomalously small coupling, since it tends to be that 
of a long cluster or bond, will imply other shorter-than- 
typical segments with larger-than-typical couplings. 

Fortunately, a modification of the naive guess of inde- 
pendence turns out to be correct, but it involves explicitly 
keeping track of the length of each segment — including 
the first and last, Ip and Z^, that link ai and ctl to the 
first and last remaining clusters,at scale F. One can show 
that a distribution of the formO 

dProh[N;AF,lF;f3i,lci;Ci,lBi ■ ■ ■ 

Cn~i,Ibn-i', PnJcni^lJl I -^,F] 

= aNiL, T)C{KF,lF)dKFdlFR{l3i,lci)d(3idlci 
P{C,i,lBi)dCidlBi X ... X C{AL,lL)dALdlL 
5{If + Ici + Ibi + ■■■ + Icn + Il-L) (25) 

for > 1 has its form preserved under renormalization 
if 



aN{L,T) = a(L,F) 



is independent of N with 



1 da 

a dT 



d/[P(C = 0,Z) + i?(/3 = 0,/)]. 



(26) 



(27) 



Note that the functions £, P and R and all other func- 
tions that will occur later are implicitly functions of F 
unless otherwise stated. 

Initially, N = L is the only term with non-zero proba- 
bility, and all the {Isi} and {Ici} are equal to ^, (a con- 
vention chosen for convenience). Therefore the needed 
solution has exactly the form of Eq (p5|). In terms of 
functions £, P, and R of the Laplace transform variable, 
y, the constraint can simply be enforced via the inverse 
Laplace transform 



dy e 



(28) 



at the end of a computation. Normalization requires that 
initially 



a(L,F,) = — 



(29) 



(which is unity in our case) independent of L. /^From 
Eqs. (F2.16) and we then have 



a(F) = - 
nr 



(30) 



for all F where np is the number density per unit length 
of the clusters remaining at scale F in an infinite system. 

The information we are interested in can be obtained 
from the function 



Z(A, L, T)dA = dProb [no cluster remaining (31) 
and - In C/L = A I L,F]. 

When the last cluster is decimated at F^, this will yield 
a contribution to Z[A = Af{Te) + Al{Te), L,Te]. In 
the Lap lace transform variable y for the lengths, using 
Eq dH) we have 



J(A,2/,F) 



LT 

L^y 



1 dZ{A,L,T) 



dT 



(32) 



= a(F) J dAF J dALS{A-AF-AL) 

R{p^O,y)C{AF,y)C{AL,y) 

since the probability of decimation of the last cluster with 
length Ic as F is increased by dV is R{(3 = 0, lc)dT. The 
function 

J{A,L,r = G)= LT-^ J{A,y,r = G) (33) 
y^L 

is then the joint probability distribution of A = — In Ci.l 
and G = - In A^;. 

The only unknown function in this process is C{A, I; F). 
It is just the joint distribution of the log-correlations of 
the first spin in a semi-infinite chain with the leftmost 
undecimated spin cluster at scale F and the length of the 
bond connecting these. The Laplace transform of C in 
L ^ y and A ^ z is simple to analyze given the additive 
form of the recursion relation Eq (pi|). This function 
C{z, y — 0) (i.e. ignoring the I variable) is exactly the 
function C{z) in F Section VB. Following the analysis 
presented there we have 



dC{z,y) 

or 



= £{z, y)[P{z, y)R{/3 = 0, y) - R{(3 ^0,y^ 0)] 

(34) 



which can be integrated. Using the special family of exact 
solutions for P and R, we obtain 



^z,y) = 



u{0)[z + ui{y)] 
ui{0)[z + u{y)] 



(35) 



with 



and 



u{y) = u{y,T), 



ui{y) = u{y,ri) 



(36) 



given generally by Eqs. (F2.51) and (F2.52) which we will 
not need in detail here. From Eqs. (F2.46) and (F2.53) 
we have 



R{f3 = Q,y)^T{y), 
defined below, and hence from Eq (|3^), 



(37) 
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J[z,y,T) = —T{y) 
nr 



z + ui[y) 



z + u{y) 

In the scaling limit of large F, small y, 
u(y) = A(y)cotli[rA(y)]-5 



(38) 



(39) 



ny) 



sinh[rA(2/)] 



where 



A(2/) = y^y + ^ 



(40) 



with lengths measured in units of ly. Since the scales 
involved initially are small, for small y we have 



ui{y) ~ w/(0). 



In the scaling limit, 



sinh'(r(5) 



(41) 



(42) 



Since in the scaling limit we expect A to typically be 
large, its distribution will be determined by the small z 
behavior, so we can also write z + ui{y) « uj{Q); [We will 
consider later what has been left out by this approxima- 
tion.] We now see that the m/(0)'s cancel in the scaling 
limit yielding 



Jiz,y,T) 



Ae 



-sr 



S + A coth(r A)] ^ sinh(r A) 



with A(?;) = \/y + 5^ and 

dProb(A, G I i) = dAdG LT-^ LT-^ J{z, y, T 

y ^ L z —f A 



our principal result. 



(43) 

G) 
(44) 



C. Distribution of Correlations 

The distribution of the log of the end-to-end corre- 
lations, A = — InGi^L, can be obtained by integrating 
Eq (jij) over G. The inverse Laplace transform in z can 
be done trivially yielding 



J(A,j/,r) = Ae*'^-^^=°*''(I^^'- 



A 



e-'^^e(A). (45) 



sinh(rA) 
At criticality, (5 = 0, we can write 

1 



sinh(rA) 



(46) 



and convert the V integral to a w integral. The inverse 
Laplace transform in y can then be done straightfor- 
wardly yielding a simple integral over w and thereby 



dProb(A \L,5 = 0) = dA Ae-AV4L 



which is of scaling form in 



A 



(47) 



(48) 



Off- critical, we have not managed to compute the dis- 
tribution of A in closed form. But progress can be made 
on the moments by expanding J in powers of z. 

In the disordered phase, for small fixed 5 > Q and small 
y, we can expand A « j/(5 + + • • •. For large L, one 
then obtains 



AP « {26L)P 



(49) 



so that the distribution is strongly peaked around 25L. 
For p = 1, Eq ( ^ ) is the first term in an expansion in 
powers of ^ / L with the correlation length 



To next order we obtain 



A: 



1 



4i 



In 



T 



+ 7 



O 



(50) 



(51) 



with 7 « 0.577, Eulcr's constant. The variance of A is 
var(A) = (A- A) 



2^ + o(L) 

£71 



(52) 



with the correction probably a power of InL. In 
Eqs. (p0|)-(p2|) we have reinserted factors of the basic 
length ly. 

Note, from Eq. (|l|) , that A is of order one for L — 
where 



26 



(53) 



is the typical correlation length. However, the true cor- 
relation length, ^ in Eq. (Eoh, is the scale at which 



A(L = 0- v/var[A(L = 0], 



(54) 



which is much larger. Thus the typical correlations are 
very small at distances of order the correlation length. 

In the ordered phase, 5 < 0, for long chains the end 
spins will almost have a spontaneous magnetization. We 
thus expect that Ci^l will be approximately the prod- 
uct of two independent end-point magnetizations of an 
infinite system. From F section VB, we thereby obtain 



ciProb[A I L] « AS^Ae-^^^^^dA. 



(55) 



This can also be computed from J using the fact that for 
6 < 0, C{z,y = 0,r ^ 00) tends to a limit as seen in F 
section VB. 
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D. Energy gaps 

We now turn to the distribution of energy gaps, more 
precisely of G = — In AE. This can be obtained from 
Eq (^) by setting z ~ 0, T = G, and inverse Laplace 
transforming over y. For general 6, we have not been 
able to perform this inverse transform analytically. But 
at the critical point the behavior again simplifies. By 
changing to the variable A ~ y/y, the inverse transform 
can be done by deforming the contour to circle a series 
of double poles, yielding, in terms of the scaled variable 



Noting that in this limit 



9 



G 

7l 



(56) 



dProb[g I L, (5 = 0] EE dgD{g, 5 = 0) 

= dg% ± e-("+i)^-V.^ + 1) 



n— — oo 
oo 



9 1 

dg^ E e-^''("+^)'(-l)™(m+-) (57) 



the second expression being the Poisson resummed form 
of the first. We see that for large g, 



while for small 

D{g,5 = 0) 



27r 



9 



(58) 



(59) 



Thus the probability of both very large and very small 
gaps is exponentially small. The former implies that the 
average gap will not decay as a power of L. Neverthe- 
less, from Eq(|5^) one can see that the integration over 
g needed to obtain the average gap is dominated by the 
maximum of the integrand yielding 



AiJ L 6 exp 




(60) 



where we have dropped a (non-universal) prefactor. 
Thus at the critical point the average gap decays 
more slowly with L — indeed with a different form of L 
dependence — than does a typical gap which has the form 
~ exp(const.L^/^). 

In the disordered phase, the behavior for L 3> ^ can be 
obtained by expanding coth(rA) in y yielding 



J{z = 0,y,T = G) 
valid for 



2de 



-25G 



[^ + 25e-^os^ 



1<^SG<^L5^ = 4- 



(61) 



(62) 



the inverse Laplace transform simply yields 
dProb[G \L,S>0] = J{z = 0,L,r = G)dG 
w 2SnGLe-^'"^dG. 



(63) 
(64) 



This result has a simple interpretationo. 

For energy scales outside of the critical region, i.e. 
r 3> the remaining clusters are typically of length 
small compared to their separations, the latter keing typ- 
ically of order l/nc- These can be thought ofS as inde- 
pendent localized excitations with density no. Thus the 
probability of one of those in a sample having an anoma- 
lously small G is proportional to the sample size L. For 
fixed L and 



(5L > G > i In ^-^ 



(65) 



Eq ( |64|) is seen to have exactly this form. Since, from 
Eq. p|), no ~ {AE)'^\ the distribution of the energies 
of a single cluster has the form 



with 



PTOh[{AE) cluster 

1 



1/z 



(66) 



(67) 



This interpretationa assumes that clusters are separated 
by of order ^ as they indeed are at the crossover scale 
away from criticality, F ~ 1/5. In this regime, length L 
corresponds to (energy) so that z plays the role of a 
continuously variable dynamic exponent as discussed in 
F section IVA and Ref. |[ Indeed, Eq (^ has exactly 
the form of a scaling function of 



(68) 



valid for j — !■ cxd and e — > with e — > any fixed value. 
Note that for extremely small gaps, G ^ SL, (for which 
the above results do not hold) the distribution depends 
on extremely rare clusters and can become less universal. 

In the ordered phase, i5 < 0, the distribution of gaps 
for L ^ ^ = ^ can be found simply. If we guess that 
G will typically be of order L, then the integral over y 
can be done by steepest descents with the saddle point 
at A{y) w jj-. The distribution of G is then found to be 
asymptotically Gaussian for L ^ 4^: 



rfProb(G \L,5 <0) 



dG 



-(G-2\5\Lf /AL 



We thus have 



and 



G = 2 |,5| L/ly 



(69) 



(70) 



varG w 2L/l^ 
where we have reinserted factors of 1^. 
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E. Relation between correlations, gaps and couplings 

We have seen that in the ordered phase the logarithm 
of the gap, — G, and in the disordered phase the logarithm 
of the end-to-end correlations, —A, both have approxi- 
mately Gaussian statistics in long samples; i.e. L ^. 
Indeed, for (5 < in this limit 



G _ 2jS\ 
L ~ 



and 



var(G) 



while similarly for 6 > 0. 

A _ 2(5 

L ly 

var(A) 2 



In J — In h 



var(ln J) + var(ln /i). 



In h — lnJ 



(71) 



(72) 



L 



ly 



= var(ln J) -|- var(ln h). 



Where does this simple behavior come from? Examina- 
tion of the basic recursion relations Eqs. ([ill) and (|l^ ) 
imply that at any scale, the effective field on any cluster 
is simply 



and the strength of any effective bond. 



J: 



with 



(73) 



(74) 



(75) 



ie section 



the sum running over all the original bonds or spins in 
the cluster (Sc) or effective bond (Ss). The last cluster 
to be decimated in a chain with free boundary conditions 
thus determines both the gap and the end-to-end corre- 
lations. If T,Ec is the sum of Eq ( [tsI) over the last to 
be decimated cluster and St is the sum over the whole 
chain, i.e. 



L-l 



St = In /j-i - In J,, 



then 



i=l 



(76) 



(77) 



since the gap is the effective field on the last cluster. 
Likewise, the log-correlations are given by 



A 



-'EC 



(78) 



since, from Eq (pT|), Ap{Te) will be the sum of In /i— In Ji 
over the section to the left of the last cluster and likewise 



AF(r_B) to the right. Eqs. ( 
the simple result that 

A- G; 



|77|) and ([78|) immediately give 



(79) 



Note that the RG procedure guarantees the physical re- 
quirements that A > and 



G > -lnf7/. 



(80) 



The specific set of bonds and spins which make up the 
last cluster for a given sample are thus very special. 

In the disordered phase for F >• | and L ^ the 
remaining clusters are typically small [of length ^ times 
factors like ln(L/^)]. Thus as L/£^ oo, 



and 



A « St + o{L) 



var At ~ var St -I- o{L). 



(81) 



(82) 



Similarly, in the ordered phase most of the system will 
be in one cluster; this is the reason for spontaneous mag- 
netization in an infinite system. Thus S ec ~ so that 



G « St + o(L) 
var G ~ var St -I- o(L) . 



(83) 



The results Eqs. (0)-(|7|) immediately follow. 

A check on this picture can be obtained by computing 
the distribution of A — G. From Eq. ( ^3|) this is found to 
be purely Gaussian with the mean and variance identical 
(within the RG approximation) to that of St, as also 
follows directly from (|7^). /^From Eqs. (76) and 
the central limit theorem tells us that the distribution of 
A — G is Gaussian with mean and variance given by 



A - G = L{A,-, 



A,,) = ?I 



var(A — G) = L(var h + var J) = — L, 



(84) 



where we used Eqs. (^) and (||) 

We conclude this section with a note on connection to 
earlier work. Shankar and Murthycl studied the correla- 
tions of quantities, such as a^af^^ or af which are local 
in fermion variables. Their results imply that, e.g.. 



In h — In J 



(85) 



with probability one for separation | r |— > cx3. Thus the 
asymptotic decay of "typical" correlations in the disor- 
dered phase of both cr^ and ct^ , as well as the local energy 
density etc., is exponential with the same "typical corre- 
lation length", ^ given by Eq. (^3|). This prediction for 
the "typical correlation length" was confirmed numeri- 
cally in Ref. ||. 
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F. Average end-to-end correlations 



The form of the average end-to-end correlation func- 
tion, Ci i, can be obtained much more simply than the 
distribution. Crudely, the average is dominated by rare 
samples in which neither of the two end spins is deci- 
mated until the final scale Te when the last cluster, which 
in this case contains both cti and ctl, is decimated. ^From 
the analysis we have carried out above, this corresponds 
simply to the 5{A) part of J(A, L, Fb). Thus up to non- 
universal factors arising from the high energy small scale 
physics, we have 

/>oo 

/ dFJ(z = 0,y,F) (86) 

^ L^yJn sinh(FA) u2(o)- 

At criticality, the integral over F is logarithmically diver- 
gent for small F yielding a ln(l/y) singularity for small 
y corresponding to 

Ci,L « § (87) 

for large L, where K \s a, non-universal constant. The 
sources of this non-universality are discussed in the next 
section. 




10 102 
L 



FIG. 2. A log-log plot of the average of the end-to-end 
correlation function at the critical point. The solid line is the 
best fit and has a slope of 0.969, close to the prediction of — 1, 
see Eq. (|^). 



IV. NUMERICAL RESULTS 



The nuigecical technique used here has been discussed 
previouslyB'BQ. It involves mapping the spin operators 
to fermions and diagonalizing therjjesulting free fermion 
problem. We refer to earlier worka-El for details. We are 
able to study sizes up to L = 256 and obtain results for 
about 50,000 samples. This large number is necessary in 
order to map out the distribution with some precision. 

A. End-to-end Correlations 

Our data for the average end-to-end correlation func- 
tion is shown in Fig. |^. The slope is about —0.97, close 
to the prediction of —1 in Eq. (p7|), and the amplitude, 
K, is close to unity. 

The distribution of values of the end-to-end correlation 
function at the critical point is shown in Fig. ^. Clearly 
the distribution is extremely broad and becomes broader 
with increasing size, even on a logarithmic scale. The 
scaling plot in Fig. ^ shows that the data fits quite well 
the analytical result of Eq. (^^. Note that there are no 
adjustable parameters in this comparison. 

Fig. 1^ shows the region near A = with a linear vertical 
scale. Whereas the data for smaller sizes appear to have 
a finite intercept at A = this decreases for larger sizes 
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FIG. 3. The distribution of the end-to-end correlation 
function at the critical point for different sizes. 
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FIG. 4. A scaling plot of the distribution of the end-to-end 
correlation function at the critical point. The solid line is the 
analytical result in Eq. (||), P(A) = (A/2) exp(-AV4). 




FIG. 5. A blowup of the region of Fig. ^ near the origin 
with a linear vertical scale. 



consistent with the analytic prediction, Eq. (|47|), of a 
linear variation for small A in the scaling limit. 

The average value in Eq. ( ^ ) arises from rare samples 
with exceptionally large values of the correlation function 
which correspond to the small A region of Figs. ^ and ||. 
We shall see that the amplitude, of the 1/L behavior 
in Eq. ( |87| ) involves both a contribution from the scaling 
function and a non-universal part which is not contained 
in the scaling function. Noting that C\,l = exp(— Ai^/^), 
the contribution to the average from the scaling function 
is 



j"^^ ^jscaling part 



1 



3-^^'''Ae-^'/4dA. (88) 



The dominant contribution is for A small, so we can re- 
place 1^ by unity, which immediately gives 



1,L 



[scaling part 



1 

2Z' 



(89) 



corresponding to an amplitude, in Eq. ( p7| ) of 1/2, 
in disagreement with the data in Fig. |^ which gives a 
value of about unity. In fact, the amplitude must clearly 
be non-universal, because one could redefine cr^ to be 
a multiple of itself which, from the scaling form in the 
logarithmic variable A, would change this amplitude but 
otherwise not change the physics. Clearly then, one 
source of the non-universal amplitude comes from the 
fact that In Ci^L in the scaling function could be replaced 
by ln(Ci^L/C'o), where Co is a constant. The scaling 
regime involves the limit InCi^^ — > — oo and so the fac- 
tor of Co is a correction to scaling (actually simply the 
leading correction to scaling discussed in F). Note that in- 
cluding the factor of Co shifts the scaling curves in Figs. ^ 
and 1^ horizontally by an amount of order 1/L^/^, thus 
giving a finite intercept at A = 0, which is seen in Fig. ||. 
Although the factor of Co is not part of the scaling func- 
tion, it does affect the average correlation function since, 
repeating the calculation which led to Eq. ( ^9|) with Co 
included, obviously gives 



1 including amplitude 
[L-l.iJav = 



Cl 

1L' 



(90) 



The factor of Co is not the only non-scaling contribu- 
tion to the amplitude. In particular, we can go back and 
look at the non-scaling parts of the distribution that arise 
from J(2, y, F) in Eq. (38). The scaling limit involved ne- 
glecting z with respect to uii/ij) in the numerator of the 
factor in parenthesis in Eq. (p8[). But for large which 
determines the behavior for small A, i.e. large Ci_l, we 
should take the quantity in parentheses equal to unity. As 
discussed earlier this corresponds to situations in which 
neither the spin cluster containing a\ nor that containing 
(Tl is ever decimated: i.e. both end spins remain active 
until the scale F^; at which the last cluster, which in this 
case contains both <Ji and cr^, is decimated. The prob- 
ability of this occurring is of order 1/i arising from the 
M/(0)'s in Eq. (^ but with a non-universal coefficient. 
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FIG. 6. A scaling plot ol the distribution of energy gaps 
at the critical point. The solid line, which fits well, is a plot 
of Eq. @ of the text. 

Thus we see that the mean correlations have various 
contributions of order 1/L each with coefficients which 
depend on the high energy small scale details that are 
not treated correctly in our RG analysis. 
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FIG. 7. Data for the distribution of ln(A£;) - ln(Ci,i), 
both at and away from the critical point, {5 = 0), showing 
that it is a Gaussian with mean and variance given by Eq. ( |8l| ) 
as expected. Note that Iv ~ 1 for the distributions used here, 
see Eq. (S). 



tions to the asymptotic form obtained from the RG equa- 
tions. 



B. Gaps 

In addition to calculating the end-to-end correlation 
function, we have also determined the energy gap at the 
critical point and show a scaling plot of the results in 
Fig. 1^. The agreement with the prediction in Eq. (|5^ ) is 
good. 

C. Results Away from Criticality 

Finally, we have also obtained some results away from 
the critical point. A particularly simple prediction is 
that ln(Ai?) — InCi.L should equal Et = Si^il^^i ~ 
^fSi^i^Ji), defined in Eq. ([79|), for each sample and 
hence its distribution should be a Gaussian with mean 
2SL/lv and variance 2L/Z„ for a// values of 6, see Eqs. ( |7^ ) 
and (0). This works very well as shown in Fig. 0, where 
the solid line is the predicted Gaussian. Indeed, we find 
that, as expected, Eq. ( [76|) holds for each sample. In 
Figs. p[]lO| we have plotted ln(Ai5) — InCi^^ against Y^t 
for L = 32 and three values of d. As expected the data 
all lie close to the curve y = x, indicated by the dashed 
line. The small deviations presumably come from correc- 



V. CONCLUSIONS 

In this paper we have studied the probability distri- 
butions of two simple properties of finite length random 
transverse field Ising chains at and near the critical point. 
In the critical regime, the natural scaling variables are 
the logarithms of the gap and the end-to-end correlation 
function divided by the square root of the length, L. The 
distributions of these properties have been computed ex- 
plicitly in the scaling limit and the results aerge well with 
numerical computations. As noted earlier^B, the aver- 
age correlation function is dominated by rare samples 
for which the end spins have correlations of order unity; 
interestingly, the information about these rare events is 
contained, up to non-universal prefactors, in the scaling 
fimction of the distribution of InCi^i/VL- 

Although the distribution of the bulk correlation func- 
tions are much harder to compute analytically, we expect 
that the same general behavior will obtain for them: At 
the critical point, A = — InCo^r/V?" is the natural scaling 
variable and the mean correlations are dominated by a 
power law singularity in the distribution at small A which 
causes the average correlations to decay as 7--(3-^)/2, 
Off critical, the typical and average correlations both de- 
cay exponentially but with the (true) correlation length 
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FIG. 8. Results for ln(A£;) - lnC(l,L) against 

St = X/i^^i ^» ~ '-^^ ^'^^ L = 32 at the critical point. 

Each point is the result for a single sample and data for 1000 
samples are shown. 



m 1 1 1 1 1 1 1 1 1 1 1 1 M 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



T 



5 = 0.203 

L=32 




30 ^ 

20 V- 
10 

^ 

-10 h 

-20 - 



-30 ^/ 

itiilllll 1 I I Illllllllllll II I w 



-30-2CM0 10 20 30 

FIG. 9. As for Fig. | but for 5 = 0.203 in the paramagnetic 
phase. 
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FIG. 10. As for Fig. | but for S = -0.178 in the ferro- 
magnetic phase. 

of the average correlations divei-ging more rapidly than 
the "typical correlation length" qH. 

This behavior is in striking contrast to that at con- 
ventional random critical points. If randomness is rele- 
vant at a classical random critical point, then the nat- 
ural variable for tjie distribution of correlations is a = 
— lnCo,r/lnr witho 



dProb(a) 



(91) 



a type of distribution sometimes called " multifractal" , 
although in terms of the log-correlations variable it 
is quite simple. The typical log-correlations are thus 
cttyp In r ± ©(Vln r) with f{atyp) = the minimum value 
of / ; i.e. they are not very broadly distributed. Nev- 
ertheless, the average correlations are dominated by rare 
pairs of spins and decay as a power law with exponent 
aav < cttyp- Off critical, the typical and average cor- 
relations decay with different correlation lengths but the 
latter is only longer than the former by a constant factor. 

The differences between conventional non-random 
quantum critical points and the random transverse field 
Ising chain are even more dramatic if one considers the 
decay of correlations in the imaginary time direction. In- 
stead of the gap scaling as a power of the length as at 
conventional quantum critical points, here the logarithm 
of the gap scales as a power of the length but with the 
ratio being broadly distributed. This implies that even 
the log-correlations in imaginary time will be extremely 
broadly distributed since they will typically behave as 
—Tin AEiocai with AEiocai loosely defined as the local 
gap. A numerical study of time-dependent correlations 
has been given in Ref. |6[ It would also be interesting to 
explore these temporal correlations in detail analytically. 
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At this point, it is natural to ask how broad is the class 
of random quantum phase transitions for which the be- 
havior is qualitatively similar to that found here. In one 
dimension, the phenomena are very general. Jadeed, one- 
dimensional random quantum Potts modclalj are in ex- 
actly the same universality class for all q and some — but 
not all — other transitions in random quantum magnetic 
chains have many similar featureaJ. 

In higher dimensions, the situation is less clear. Sim- 
ulations o n I . qu antum spin glasses in two and three 
dimensionsElrlia have suggested fairly|-eoB|Ventional scal- 
ing. However, recent simulationaid^a of the two- 
dimensional random quantum Ising ferromagnet, for 
which larger sizes and better statistics can be obtained, 
indicate that most of the features found in one-dimension 
also apply to this case. The possibility that some strongly 
random quantum critical points in higher dimensions 
share some of the exotic features of the random trans- 
verse field Ising chains is quite intriguing. 
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APPENDIX A: 

In this appendix, we briefly consider the distribution of 
the gap in long chains with periodic boundary conditions 
like those used in Ref. |[ Note that in such systems, as in 
bulk systems, the distributions of correlations are much 
more difficult to compute analytically than the end-to- 
end correlations for the free boundary case and we have 
no new results on these. 

With periodic boundary conditions, it can be shown 
that, given N remaining clusters at scale F in a system of 
length L, the joint distribution of the effective couplings 
and their lengths is simply the conditional distribution 
of N otherwise independent clusters and bonds in an in- 
finite system given that their total length is L. As for 
the free boundary case, this implies correlations between 
the couplings which can be handled by Laplace transform 
methods as in the main text. The probability that there 
are N clusters at scale F can be computed from the RG 
equation that couples the flow of the A^-cluster distribu- 
tion and that of — 1 clusters when one of the N bonds 
(or clusters) is decimated. As in the free case, G, minus 
the logarithm of the gap, is determined by the scale at 
which the last cluster is decimated. We obtain 



dProb(G) dG LT-^ |_ 
Fj — > II si 



L ^ y 1^ sinh^(GA) [A coth(GA) - 6] 

G ^2 

dV- 



-(G-r)[A coth(rA) + 







sinh^(FA) 



(Al) 



We have not been able to perform the integrations in 
Eq. (Al) analytically, but the asymptotic behavior of the 
scaling distribution at the critical point can be found. In 
terms of the rescaled variable g — G/^/Tj, we find for 
small g 



dFvoh{g) « dg 



2g3 



very similar to the free-end case Eq. (| 
dProb(g) 



do —e 
^ 2 



(A2) 

For large 5, 
(A3) 



again quite similar to the free-end result Eq. (p9|). These 
limiting forms are a reasonable fit to the numerical data 
of Rcf. 0. 
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